TeV Neutrinos from Companion Stars of Rapid-Rotating Neutron Stars 
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We have estimated the number flux of of 's which are produced due to the hadronic interactions 
between the cosmic rays coming from a neutron star and the matter in a companion star. The event 
rate at 1 km 2 detectors of high-energy neutrinos such as ICECUBE, ANTARES and NESTOR is 
also estimated to be 2.7 x 10 4 events yr _1 when the source is located at 10 kpc away from the 
Earth. We have estimated the number of such a system and concluded that there will be several 
candidates in our galaxy. Taking these results into consideration, this scenario is promising and will 
be confirmed by the future observations. 
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Neutron stars have been suggested as possible sources 
of very high energy cosmic rays, including the ultra- 
high-energy cosmic rays whose energies are greater than 
10 20 eV (see Refs. jl| and references therein). It is gener- 
ally considered that a magnetosphere is formed around a 
neutron star where the number density has the Goldreich- 
Julian value II. It is also pointed out that the condi- 
tions which lead to the conversion of most Poynting flux 
into kinetic energy flux in cold, relativistic hydromag- 
netic winds exist around the magnetosphere (i.e. the 
wind zone) ||. 

On the other hand, it is considered that most of the 
normal stars are born in binaries. Observationally, it is 
reported that somewhere between 50 and 100% of the 
normal stars are in binaries [@. Thus it is very natu- 
ral to consider the situation that a binary is composed 
of a neutron star which emits high energy cosmic rays 
and a red super giant whose radius is ~ 10 14 cm. What 
happens if the high energy cosmic rays rush into the com- 
panion? This is the theme in this paper. It is noted that 
the wind from the supergiant is much weaker than the 
pulsar wind and should be blown off. Thus it is not nec- 
essary to consider the effects of interaction between the 
wind from the companion star and the pulsar wind. 

The conclusion is as follows. Charged pions are pro- 
duced due to the hadronic interactions between the cos- 
mic rays and the stellar matter. Before losing their en- 
ergies, the pions can decay into leptons, including high 
energy (TeV) v^'s. The v^'s can be detected by 1 km 2 
detectors of high-energy neutrinos such as ICECUBE, 
ANTARES and NESTOR §. The expected event rate 
is 2.7 x 10 4 events yr _1 . Also, we have estimated the 
number of such a system and concluded that there will 
be several candidates in our galaxy. It is noted that high 
energy neutrinos from Cygnus X-3 were calculated on the 
basis of the binary model in some excellent papers Q , but 
the companion star was assumed to be a main-sequence 
star whose column density seems to be too high to pro- 
duce high energy neutrinos (see Eq. |l5|). In this study, 
the companion star is assumed to be a red super giant so 
as to solve the difficulty and neutrino flux is estimated 
using the result of the theoretical accerelation model of 
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cosmic rays around the magnetosphere 0, @- Also, we 
estimate the number of candidates in our galaxy. 

From now on, the energy spectrum of high energy cos- 
mic rays is reconsidered. Then the energy loss processes 
for the cosmic rays and their daughters in a companion 
star are investigated. Next, the number flux of v^s at 
the Earth from a single source is estimated. Finally, the 
event rate and number of candidates are discussed. 

Here the energy spectrum of high energy cosmic 
rays is estimated (see Ref. |Q, for details). The 
Goldreich- Julian value can be expressed as nQj(r) = 
B(r)£l/ (AirZec), where B(r) is the strength of the mag- 
netic field, is the angular velocity of the neutron star. 
It is assumed that the magnetic field in the magneto- 
sphere is dominated by the dipole component and can 
be expressed as B(r) = B s (R s /r) 3 , where B s and R s 
are the strength of the surface magnetic field and radius 
of the neutron star, respectively. In this study, we set 
R s = 10 km throughout. The light cylinder radius can 
be expressed as R\ c = 10 7( 1^ cm, where tt 3 k is Sl/3000 
rad s — 1 . The typical energy of the accelerated cosmic 
rays E cr can be estimated by assuming that the equipar- 
tition of energy is realized between the magnetic fields 
and cosmic rays (E'er = -B(r) 2 /87TOGj(r)). At R\ c in the 
equatorial plane, E cv can be written as 



E clfi = 1.5 x 10 9 ZB sA2 nl k GeV 



(1) 



where i? s ,i2 is -B s /10 12 Gauss. The strength magnetic 
field at (i?i c , h) in the cylindrical coordinate, where h 
is the height measured from the equatorial plane, can 
be expressed as B(R\ c ,h) ~ B{{R\ C , 0)) cos a, where 
cos a — Ric/ (Rf c + ft 2 ) 1 / 2 . Thus E CI depends on h and it 
can be expressed as E cr (h) ~ -E C r,o c °s 3 ex. Thus the en- 
ergy spectrum of high energy cosmic rays from the light 
cylinder can be expressed as 
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dE * dt Z{E CTfi E*y/z [(E CI , Q /E z y/* - 1] 
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(2) 



where E z represents the energy of the cosmic rays whose 
charge number is Z . In this study, we consider two cases 
for the chemical composition of the cosmic rays, that is, 
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protons and irons. The chemical composition of the cos- 
mic rays may be dominated by protons, which will be 
ejected from the companion star and the magnetosphere 
may be filled with the protons. On the other hand, the 
chemical composition of the cosmic rays may be domi- 
nated by irons, which will be ejected from the surface of 
the neutron star. 

Since the energy spectrum of cosmic rays can be ob- 
tained, we consider the energy loss processes for the cos- 
mic rays and their daughters in a companion star. At 
first, we assume that the chemical composition of the 
cosmic rays is dominated by protons. We consider the en- 
ergy losses due to ionization, synchrotron radiation, pair 
production due to proton-electron interaction, photopion 
production, proton-proton interaction, respectively. In 
order to calculate the energy loss timescale, the condi- 
tions in the companion star has to be specified. Here we 
assume that the radius and temperature of the red su- 
per giant is 10 14 cm and 3000K, respectively. Thus the 
number of photons in the star can be estimated as 



n 7 = 5.4 x 10 11 



(— V 

V3000K7 



(3) 



The average baryon number density is used in this 
study for simplicity and is estimated as 



n B = 2.8 x 10 



is f M 



V 10M e 



10 14 cm 



-3 



(4) 



where i? s t a r is the radius of the companion, M is the mass 
of the companion and M & is the solar mass. 

Now we investigate the energy loss timescale. As for 
the energy loss rate due to coulomb interactions, colli- 
sions are dominated by the thermal electrons ||. Thus 
the energy loss rate can be written approximately as 



dE\ 30ca T m e c 2 Z 2 



dt 
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n e eV s , 



(5) 



where ctt is the Thomson cross section, (3 is the velocity 
of the cosmic ray divided by the light speed. Thus the 
energy loss timescale can be estimated as 



3.2 x 10 b 
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which is, as we show below, too large to be effective. 

As for the synchrotron energy loss timescale, it can be 
written as 



"p,syn 



X te,syn ^ IT x 10 t e S y n , (7) 



where i e ,syn is the synchrotron energy loss timescale 
for an electron and can be written as t etSyn = 3.9 x 

10 (1 GeV/E) (10 2 G/B) 2 s. Thus synchrotron energy 
loss is negligible. 



As for the pair production process p + 7 — > p 
e + , the energy loss timescale can be written as 
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where a is the fine-structure constant, r e is the classi- 
cal electron radius, T is the lorentz factor of proton and 
temperature is assumed to be 3000K. This timescale is 
too large to be effective in the energy range where we are 
interested in. 

As for the photopion process, the A-resonance occurs 
at fl 



,,, = 2 |um1?X] G.-V. 



(10) 



At the A-resonance, the cross section becomes about 
5xl0~ 28 cm~ 2 jl(J. Even if this value is adopted, the 
interaction length for this process can be estimated as 



1 



ll5 /3000K\ 



J>7 



5 x 10" 28 n. 
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7x 10^^——j cm. (11) 



Thus energy loss timescale due to photopion production 
is 



^ 3000K y 



(12) 



On the other hand, the cross section of pion produc- 
tion due to proton-proton interaction is about 50mb [0 . 
Thus the energy loss timescale due to proton-proton in- 
teraction is 
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(TppTlBC 



2.4 x 10 



Km) \10 14 cm 



(13) 



This timescale is much shorter than the other timescales. 
Also, the interaction length l pp = ct pp is much shorter 
than the radius of the companion. 

From the above discussion, we conclude that the domi- 
nant energy loss process is proton-proton interaction and 
charged pions are produced due to the interactions. The 
high energy protons lose ~ 50% of their energies per each 
interaction. Here we have to investigate whether these 
pions can decay into leptons before losing their energy. 
If they can, the leptons also energetic and the energy of 
each neutrino is about 25% of the pion's energy. Since 
the lifetime of a charged pion in its rest frame is 2.6 x 10~ 8 
s, the observed lifetime is 
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10 15 eV 



s. 



(14) 



On the other hand, the cross section of proton-pion inter- 
action is about 30mb |1Q] . Thus the energy loss timescale 
due to this interaction can be estimated as 
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(15) 



H=3xlO-( T ^) 2 ^ eV ^2g/(^)P(^) (19) 



Thus, charged pions with energy below ~ 2 x 10 15 eV 
can decay into leptons before losing their energy. 

We also check the opacity of the companion star for 
the high energy neutrinos. The neutrino- nucleon total 
cross section can be written approximately as ct„jv ~ 5 x 
10 -39 [E v j\ GeV) cm 2 11 1. Thus the average number 



2.7 x 10 4 



1 km 



Bli 3 nl{ 3 Z- 5/3 events yr" 1 , (20) 



of interaction for a single neutrino in the companion star 
can be estimated as 



I N = 1.4 x 10" 



M 



10 14 cm\ / E v \ 



1OM J ^ i? star ) \lGeVj 



(16) 



From this estimate, we can regard that the neutrinos can 
escape from the companion freely. 

Now we can estimate the number flux of v^s [s -1 cm~ 2 
GeV -1 ] at the Earth from a single source. It is assumed 
that the source is located at 10 kpc away from the Earth 
(D = lOkpc). The expected number flux is (E v < 5 x 
10 14 eV) 



f(E v )= 
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dN v 



AttD 2 dEudt 



N 2 (E p ) 12dN p 



8.5 x 10 



_ 17 / lOkpc 
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4ttD 2 dtdEn 



R l/3 2/3 7 - 5 



where d is the diameter of the detector. The factor 2 in 
the integral is introduced to count in both of contribu- 
tions from v^s and v^s. Strictly speaking, the factor 
(5.6/47T ~ 0.45; see Eq. (|27j)) which takes into account 
the solid angle measured from the companion star should 
be multiplied in Eq. (]2(]). Since the number of detected 
atmospheric neutrino background by the planned neu- 
trino telescopes whose angular resolution 9 is about 1° 
is estimated as N hkg ~ 10 (0/ldeg) 2 km -2 yr -1 @, the 
signal from the companion star will be easily detected. 

We add a comment on the chemical composition of 
the high energy cosmic rays. In the above discussions, 
the chemical composition is assumed to be proton. In 
the case that the chemical composition is iron, we have 
to estimate the fragmentation timescale. There is an em- 
pirical formula for the total inelastic cross-section of pro- 
tons on nuclei with mass number A > 1 |15|. Accord- 
ing to the formula, the cross section can be written as 
cr tot = 45,4 a7 [l + 0.016 sin(5.3- 2.63 InA)] mb. Thus 
the fragmentation timescale can be estimated as Tf = 
5/3 N z (E p ) ^ 3 g x 1Q -2 (56/A) ' 7 (io 15 cm- 3 /n B ) s. This timescale is 

eV z shorter than the proton-proton interaction timescale (see 



(17) 



where N(E p ) = 0.61 + 0.56 ln(1.88S p ) +0.129 ln 2 (1.88£:p) 
is mean value of the charged-particle multiplicity as a 
function of ^fs JlJ] . It is noted that the mean energy of 
can be written as E P /12N(E P ) when the average number 
of 7T°, 7r + and 7r - produced per one inelastic interaction is 
same with each other. As discussed above, there may be 
a cut-off around £„~5x 10 14 eV due to the interaction 
of the charged pions and stellar matter. Of course, such 
a cut-off energy strongly depends on the radius of the 
companion (see Eq. |l^). The neutral pions produced 
by pp collision decay into gamma-rays, which may also 
produce high energy pions through the pj interactions. 
Thus the present estimation for the flux of neutrinos will 
be a lower limit. Detailed Monte Carlo calculation is 
needed to obtain the precise flux of neutrinos, which is 
now underway and is presented in the forthcoming paper. 
It is also noted that the high energy gamma-rays can not 
penetrate the companion star and lose their energy by 
interactions with the matter in the companion. 

Finally, the event rate for a single source at 1 km 2 de- 
tectors of high-energy neutrinos is estimated. The prob- 
ability that a muon neutrino will produce a high-energy 
upward moving muon with energy above 2 GeV in a ter- 
restrial detector is approximately given by the ratio of 
the muon range to the neutrino mean free path fl3]| . The 
result can be fitted as P{E V ) = 1.3 x 10 -9 (E v /l TeVf 
where /? = 2 for E v < 1 TeV and (3 = 1 for E v > 1 TeV. 
Thus the expected event rate R [yr _1 ] is (over 27T sr) 



shorter than the proton-proton interaction timescale (see 
Eq. 13). Even if the chemical composition is composed of 
deuterium (A = 2), the fragmentation timescale is com- 
parable with the proton-proton interaction timescale as 
long as riB ~ 3 x 10 15 cm~ 3 . Thus we can conclude that 
the fragmentation occurs at first when the cosmic rays 
rush into the companion star and subsequent reactions 
are same with the case in which chemical composition of 
the cosmic rays is proton. As for the number flux of nu- 
cleus after the fragmentation, it will be (A/Z) 5 ^ 3 ^3.2 
times larger than the proton case (see Eq. g). Thus the 
expected event rate of high energy neutrinos will be also 
enhanced by a factor of ~ 3.2. 

When the neutrino mixing angles are really large [ fl6| , 
the fluxes of v^'s and v^s will fall off by a factor of 3. 
In this case, ^ T 's might be detected as the double bang 
events Jl7j ]. 

We also present a rough estimate for the number of 
candidates (N) in our galaxy. It may be expressed in the 
form 



A = r x SNR x / x I — 
1 4tt 



(21) 



where r is the duration [yr] for a source to emit the 
high-energy f M 's, SNR is the supernova rate [yr -1 ] in 
our galaxy at present, / is the probability that a mas- 
sive star which forms a neutron star at its death has a 
companion whose mass is above ~ 8M Q , dfl is the aver- 
age solid angle of a companion star measured from the 
neutron star. 
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The duration r will be given as r = min(r sta r, Tns)j 
where r s t ar is the lifetime of the massive companion star 
(~ 10 6 yr) and tns is the duration in which the neutron 
star can emit the high energy cosmic rays. From Eq. |TJ, 
the energy of the cosmic rays should be greater than ~ 
10 eV in order to produce much energetic neutrinos. 
On the other hand, the angular velocity of a neutron 
star decreases as 



i 



rad s 



(22) 



where £li is the initial angular velocity and I ~ 10 45 g 
cm 2 is the inertial moment of a neutron star. From Eq. [l], 
E C i,o becomes ~ 10 16 eV when B Sy i2 = 1 and f^fe = 0.1. 
Thus, if f^i is much larger than 0.1, tns can be estimated 
as 



T NS ~ 3 X 10' 



/ 10 12 G 



yr, 



(23) 



which is shorter than r s t a r- Thus it is concluded that r 
= tns- 

The present SNR is about 10~ 2 yr" 1 @. As for /, 
since most of the normal stars are considered to be born 
in binaries we can estimate this value only by using 
the salpeter's initial mass function. The result can be 
written as |19fl 



/a 25 dm 



-2.35 



f,, 1 ? 5 dm m 2 35 



2.4 x 10~ 2 . 



(24) 



Finally, the average solid angle of a companion star 
measured from the neutron star is estimated. Although 
there is much uncertainty, Popova et al. (1982) discussed 
that the distribution of the distance (a) between the bi- 
nary stars may obey the relation as follows: 



dN 
da 



A 



cm 



lOi? < a < 10 b R 



(25) 



where A = (5 In 10) _1 is the normalization factor and i? 
is the solar radius. Since the radius of the companion 
is assumed to be 10 14 cm, the opening angle 8 can be 
estimated as sin# = 10 14 /a as long as a is greater than 
10 14 cm. The solid angle can be estimated as 



dQ 



2tt 



sin 9d9 = 2tt 



10 2 



28 \ 1/2' 



(26) 



In case a < 10 14 cm, we assume that the solid angle is 
An. We can easily show by Eqs. (p|)-(p^|) that the domi- 
nant energy loss process is still proton-proton interaction 
within t ~ tns even if the pulsar wind blows off the en- 
velope of the companion star. Now we can estimate the 
mean solid angle of a companion measured from the neu- 
tron star. It can be calculated as 



f 10 dN 
dtt = An—^da 

JWRq 
rW 6 R e 

2tt 

'10 14 



da 
1 



(27) 
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da ~ 5.6 



At last, we can estimate the number of candidates (N) 
in our galaxy. It is written as 

5.6 



N ~ 3 x 10 4 x KT 2 x 2.4 x 10~ 2 x 



47T 



= 3.2. 



(28) 



We have estimated the number flux of of v^s which are 
produced due to the hadronic interactions between the 
cosmic rays coming from a neutron star and the matter 
in a companion star. The event rate at 1 km 2 detectors 
of high-energy neutrinos such as ICECUBE, ANTARES 
and NESTOR is also estimated to be 2.7x 10 4 events yr" 1 
when the source is located at 10 kpc away from the Earth. 
Finally, we have estimated the number of such a system 
and concluded that there will be several candidates in 
our galaxy. Taking these results into consideration, this 
scenario is very promising and will be confirmed by the 
future observations. 
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